Effect of local pumping on random laser modes 
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We have developed a numerical method based on the transfer matrix to calculate the quasimodes 
and lasing modes in one-dimensional random systems. Depending on the relative magnitude of the 
localization length versus the system size, there are two regimes in which the quasimodes are distinct 
in spatial profile and frequency distribution. In the presence of uniform gain, the lasing modes have 
one-to-one correspondence to the quasimodes in both regimes. Local excitation may enhance the 
weight of a mode within the gain region due to local amplification, especially in a weakly scattering 
system. 
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I. INTRODUCTION 



Random laser, in which optical feedback is provided 
by scattering of light due to spatial inhomogeneity of the 
medium rather than by well denned mirrors, has recently 
attracted much attention^. One important topic of re- 
search is the nature of random laser modes. For a random 
laser with non-resonant feedback, the lasing modes are 
the diffusive modes, i.e, the eigenmodes of the diffusion 
equation 2 -. For a random laser with resonant feedback, 
the lasing modes are believed to be the quasimodes, i.e, 
the eigenmodes of the Maxwell equations 3 -. This belief 
implies the quasimodes of a passive random system are 
not modified by the presence of gain. Such assumption 
is confirmed by the numerical studies of lasing modes in 
the localization regime^. With the introduction of gain, 
the localized modes of a passive random system are pre- 
served and serve as the lasing modes. This conclusion is 
extended to the random systems far from the localization 
regime without direct confirmation. The lasing modes 
are regarded as the quasimodes with small decay rate, in 
particular the anomalously localized states^. However, 
recent theoretical study£ reveals that the quasimodes of a 
passive random system are not the genuine normal modes 
of the same system with gain. This is because the spatial 
inhomogeneity of dielectric constant introduces a linear 
coupling between the quasimodes, mediated by the po- 
larization of the active medium. The latest development 
of semiclassical laser theory for open complex or random 
media leads to the speculation that the lasing mode in 
a weakly scattering system may be a composite of many 
quasimodes with low quality factor— Moreover, under 
local excitation the reabsorption outside the local gain 
region suppresses the feedback from the unpumped part 
of the random sample and effectively reduces the sys- 
tem sizeii. The lasing modes are therefore completely 
different from the quasimodes and confined in the vicin- 
ity of the pumped region. All these studies prompt us 
to investigate carefully the relation between the lasing 
modes and the quasimodes in both global pumping and 



local pumping. In this paper, we address the question 
whether the lasing modes are the quasimodes of pas- 
sive random systems. The answer to this question de- 
termines whether the statistical distribution of the decay 
rates of quasimodes can be used to predict the lasing 
threshold and the number of lasing modes for random 
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We conduct detailed numerical studies of quasimodes 
and lasing modes in one-dimensional (ID) random sys- 
tems. A numerical method based on the transfer ma- 
trix is developed to calculate the quasimodes as well 
as the lasing modes in the presence of global or local 
gain. The main advantage of this method as compared 
to the finite-difference time-domain method is that it 
can calculate the quasimodes of weakly scattering sys- 
tems which overlap spectrally and have short lifetime. 
In our numerical simulation, the scattering strength is 
varied over a wide range. The quasimodes, as well as 
the lasing modes, are formed by distributed feedback in 
the random system. The conventional distributed feed- 
back (DFB) laser, made of periodic structures, operates 
either in the over-coupling regime or the under-coupling 
regime 2 ^. The random laser, which can be considered as 
randomly distributed feedback laser, also has these two 
regimes of operation, fn the under- coupling regime the 
system size L is much less than the localization length £, 
while in the over-coupling regime L > £. The dominant 
mechanism for the mode formation differs in these two 
regimes, leading to distinct characteristics of mode pro- 
file and frequency distribution. With the introduction 
of uniform gain, the lasing modes have one-to-one corre- 
spondence to the quasimodes in both regimes. However, 
local pumping can make the lasing modes significantly 
different from the quasimodes, especially in the under- 
coupling systems. Some quasimodes even fail to lase no 
matter how high the pumping level is. The results we 
obtain help understanding the random lasing with reso- 
nant feedback in the weakly scattering systems 2 ^, espe- 
cially the recent observations of periodic lasing peaks in 
frequenc y 22 ' 23 . 
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II. NUMERICAL METHOD 

We have developed a numerical method based on the 
transfer matrix to compute the quasimodes of ID passive 
systems. This time-independent method is also applied 
to the calculation of lasing modes at the threshold under 
global or local excitation. The random system is ID lay- 
ered structure. It is composed of N dielectric layers with 
air gaps in between. The refractive index of the dielec- 
tric layers is rid, and that of the air gaps is 1. Both the 
thickness d\ of the dielectric layers and the thickness of 
air gaps di are randomized. d\.2 = ^1,2(1 + o~v)i where 
< a < 1 represents the degree of randomness, and r\ 
is a random number in [—1,1], d\ (cfe) is the average 
thickness of the dielectric layers (air gaps). Outside the 
random system the refractive index is constant and its 
value is equal to the average refractive index n e f f of the 
random system to eliminate the boundary reflection. 

According to the transfer matrix formula: 




where po and qo represent the forward and backward 
propagating waves on one side of the random system, 
pi and qi on the other side, M is a 2 x 2 transfer matrix 
that characterizes wave propagation through the random 
system. The eigenmode of such an open system can be 
defined as "natural mode" or "quasimode" , which gener- 
alizes the concept of eigenmode of a closed system^. It 
satisfies the boundary condition that there are no incom- 
ing waves but only outgoing waves through the boundary 
of a random system, namely po — and q± = 0. In a pas- 
sive system (without gain or absorption, the refractive 
indices being real numbers), such boundary condition 
requires the vacuum wavevector be a complex number, 
kg = ko r + ikoi. Substituting the boundary condition 
into Eq. [TJ we get M22 = 0. Since M22 is a complex 
number, both the real part and imaginary part of M22 
are equal to 0. These two equations are solved to find ko r 
and feoj. kg r — uj/c tells the frequency w of a quasimode, 
and koi — —^/c gives the decay rate 7 of a quasimode. 

After finding ko of a quasimode, the correspond- 
ing wavefunction can be obtained by calculating the 
electric field distribution E(x) throughout the random 
system with the transfer matrix M(ko). The wave- 
function inside the random system can be written as 
E(x) = E+(x)e m ^ k ° x + EJ(x)e- m ^ k ° x , where n(x) 
is the (real part of) refractive index at position x, 
E+(x)e m ( x ^ k ° x represents the forward-propagating field, 
and E_{x)e~ m ( x ' k ° x the backward-propagating field. 
Since kg is a complex number, the amplitudes of forward 
and backward propagating fields are E + ( x )e~ n ( x ' koiX an j 
E_(x)e n ( x ' koiX (k 0i < 0). These expressions show that 
there are two factors determining the wavefunction. The 
first is E±(x), which originates from the interference 
of multiply-scattered waves. The second is e ±n ^ koiX , 
which leads to exponential growth of the wavefunction 



toward the system boundary. Outside the random sys- 
tem, the wavefunction grows exponentially to infinity due 
to the negative kot- This is clearly unphysical. Thus we 
disregard the wavefunction outside the random system 
and normalize the wavefunction within the random sys- 
tem to unity. 

Optical gain is introduced to the random system by 
adding an imaginary part rii (negative number) to the 
refractive index. In the case of uniform gain, n, is con- 
stant everywhere inside the system. Outside the random 
system rii is set to zero. Different from the quasimode 
of a passive system, the vacuum wavevector ko of a las- 
ing mode is a real number. The wavevector inside the 
random system is a complex number, k = k r + iki — 
ko[n(x) + irii]. Its imaginary part ki = k^rii is inversely 
proportional to the gain length l g . The onset of lasing os- 
cillation corresponds to the condition that there are only 
outgoing waves through the boundary of the random sys- 
tem. The absence of incoming waves requires M22 — in 
Eq. (1). Again since M22 is a complex number, both its 
real part and imaginary part are zero. These two equa- 
tions are solved to find ko and rij. Each set of solution 
(ko, n,) represents a lasing mode, ko = uj/c sets the lasing 
frequency u>, and riiko = ki = l/l g gives the gain length 
l g at the lasing threshold. The spatial profile of the las- 
ing mode is then obtained by calculating the field distri- 
bution throughout the random system with the transfer 
matrix M(ko, rii). Since our method is based on the time- 
independent wave equation, it holds only up to the lasing 
threshold^. In the absence of gain saturation, the ampli- 
tude of a lasing mode would grow in time without bound. 
Thus we can only get the spatially-normalized profile of a 
lasing mode at the threshold. The lasing mode is normal- 
ized in the same way as the quasimode for comparison. 
The amplitudes of forward and backward propagating 
fields of a lasing mode are E + (x)e~ nik ° x and E- (x)e nik ° x 
(rii < 0). The exponential growth factors e ±riik ° x depend 
on the gain value |njfco|- 

Local pumping is commonly used in the random laser 
experiment. To simulate such situation, we introduce 
gain to a local region of the random system. Our method 
can be used to find the lasing modes with arbitrary spa- 
tial distribution of gain. The imaginary part of the re- 
fractive index rii(x) = hif(x), where f(x) describes the 
spatial profile of gain and its maximum is set to 1, hi 
represents the gain magnitude. The lasing modes can 
be found in the way similar to the case of uniform gain. 
The solution to M22 — gives the lasing frequency ko 
and threshold gain n^fco. The normalized spatial profile 
of a lasing mode is then computed with M(k , hi). 



III. RESULTS AND DISCUSSION 

Using the method described in the previous section, 
we calculate the quasimodes of ID random systems. 
The quasimodes are formed by distributed feedback from 
the randomly-positioned dielectric layers. We investi- 
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gate many random structures with different scattering 
strengths. Depending on the relative values of the lo- 
calization length £ and the system length L, there are 
two distinct regimes in which the quasimodes are dra- 
matically different: (i) over-coupling regime L > £; (ii) 
under-coupling regime L <C £. 

As an example, we consider the random structure with 
d% = lOOnm and di = 200nm. a = 0.9 for both d\ 
and c?2- To change from the under-coupling regime to 
the over-coupling regime, we increase the refractive in- 
dex rid of the dielectric layers. In particular, we take 
rid = 1-05 and 2.0. The larger rid leads to stronger scat- 
tering and shorter localization length £. To obtain the 
value of £, we calculate the transmission T as a func- 
tion of system length L. (InT) is obtained from aver- 
aging over 10,000 configurations with the same L and 
a. When L > £, (lnT(L)) decays linearly with L, and 
= —d(hiT(L))/dL. In the wavelength (A) range of 
500nm to 750nm, £ exhibits slight variation with A due 
to the residual photonic bandgap effect. For rid = 1-05, 
£ - 200 - 240^m, while for n d = 2.0, £ - 1.2 - 1.5/im. In 
the calculation of quasimodes, we fix the number of di- 
electric layers N = 81 and (L) = 24.1/im. For n = 1.05, 
£ 3> L in the wavelength range of interest, thus the ran- 
dom system is in the under-coupling regime. In contrast, 
for n — 2.0, £ <C L and the system in the over-coupling 
regime. 

To illustrate the difference between over-coupling 
regime and under-coupling regime, we compare the quasi- 
modes of the same random structure with different n d , 
namely, nd = 2.0 or 1.05. Figures 1(a) and (b) are 
the typical transmission spectra of these two systems. 
For the system with nd — 2.0 most transmission peaks 
are narrow and well separated in frequency, while for 
nd = 1.05 the transmission peaks are typically broad and 
overlapped. We find ko = feor + ikoi of the quasimodes in 
the wavelength range of 500-750nm. Figure 1(c) shows 
the values of kg r and fcoi/(fcoi) of these modes ((feoi) is the 
average over all the quasimodes in the wavelength range 
of 500-750nm). In the system with nd = 2.0, most quasi- 
modes are well separated spectrally, and they match the 
transmission peaks, ko r corresponds to the frequency of 
a transmission peak, and k$i to the linewidth of a trans- 
mission peak. However, some quasimodes are located 
close to the system boundary, thus having relatively large 
koi. They are usually invisible in the transmission spec- 
trum due to spectral overlap with neighboring transmis- 
sion peaks, which cause the number of transmission peaks 
[Fig. 1(a)] slightly less than the number of quasimodes 
[solid squares in Fig. 1(c)]. In the system with nd — 1.05, 
however, the number of peaks or maxima in the trans- 
mission spectrum [Fig. 1(b)] is significantly less than the 
number of quasimodes [open circles in Fig. 1(c)]. This 
is because in the under-coupling regime the decay rates 
of the quasimodes often exceed the frequency spacing to 
neighboring modes. The spectral overlap of the quasi- 
modes makes the transmission peaks less evident and 
some even buried by the neighboring ones. It is clear in 
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FIG. 1: (Color online) (a), (b): Transmission T through a 
ID random structure with nd = 2.0, 1.05 as a function of vac- 
uum wavevector ko- (c) Frequencies ko r and normalized decay 
rates koi/{koi) of the quasimodes in the random systems with 
rid = 2.0 (solid square) and rid = 1-05 (open circle), com- 
pared with the frequencies ko and normalized threshold gain 
hi /(hi) of lasing modes in the same systems with nd = 2.0 
(+) and nd = 1-05 (x) under uniform excitation, (d) Nor- 
malized frequency spacing Afco r / (Afcor) of neighboring quasi- 
modes in the random systems with rid = 2.0 (solid square) 
and nd = 1-05 (open circle), compared with the normalized 
frequency spacing Afco/(Afco) of neighboring lasing modes in 
the same systems with nd = 2.0 (+) and rid = 1-05 (x) under 
uniform excitation. 



Fig. f (c) that the decay rate fluctuation is much stronger 
in the random system with n d = 2.0 (solid squares) than 
that with nd = 1.05 (open circles). This is consistent 
with the broadening of quasimode decay rate distribu- 
tion as a system approaches the localization regime with 
increasing scattering strength. Figure 1(d) plots the fre- 
quency spacing Afco r between adjacent quasimodes nor- 
malized to the average value (Afcor)- The quasimodes 
of the random system with nd — 1.05 are more regu- 
larly spaced in frequency than those in the system with 
nd = 2.0. The average mode spacing is inversely propor- 
tional to the system length L. 

To interpret this phenomenon, we investigate the wave- 
functions of the quasimodes. Figure 2(a) [(b)] shows 
the spatial distribution of intensity I(x) = [^(a;)! 2 for 
a typical quasimode of the random system with n d = 2.0 
[rid = f.05). I{x) is normalized such that the spatial in- 
tegration of I(x) within the random system is equal to 
unity. The expression of E(x) given in the previous sec- 
tion reveals the two factors determining the envelop of 
the wavefunction, i.e., the interference term E^(x) and 
the exponential growth term e ± n ( x ) k oix_ Depending on 
which term is dominant, the spatial profile of the quasi- 
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FIG. 2: (Color online) Spatial intensity distributions of quasi- 
modes (black solid curve) and the corresponding lasing modes 
in the presence of global gain (red dashed curve) or local gain 
(green dotted curve). The pumped region is between the two 
vertical lines, L p — 11.87/im. (a) The mode marked as A in 
Fig. 1(a), A = 646nm, rid = 2.0. (b) The mode marked as B 
in Fig. 1(b), A = 549nm, n d = 1.05. 



modes can be drastically different. In the over-coupling 
regime, strong scattering makes the interference term 
dominant, and I(x) exhibits strong spatial modulation. 
Most quasimodes are localized inside the random system, 
similar to the mode in Fig. 2(a). Their decay rates are 
low as a result of the interference-induced localization. 
In the under-coupling regime, the interference effect is 
weak due to small amount of scattering. The exponen- 
tial growth term e ±"( a ) fe oix d om ma,tes E(x), making I(x) 
increase exponentially towards the boundaries. The in- 
terference term only causes weak and irregular intensity 
modulation. A typical example of such mode profile is ex- 
hibited in Fig. 2(b). Since the quasimodes in the under- 
coupling system are spatially extended across the entire 
random system, the rates of light leakage through the 
boundaries are much higher than those of the localized 
modes in the over-coupling system. We repeat the above 
calculations with many random systems, and find the two 
different types of quasimodes are rather typical for the 
systems in the over-coupling and under-coupling regimes. 
The mode profiles and frequency spacings in the under- 
coupling systems reveal the feedbacks from the dielectric 
layers close to the boundaries are dominant over those 
from the interior. Thus the quasimodes in the under- 
coupling systems are formed mainly by the feedbacks 



from the scatterers near the system boundaries. How- 
ever, the feedbacks from the scatterers in the interior 
of the system are weak but not negligible, e.g. they in- 
duce small fluctuations in the frequency spacings and the 
decay rates. Note that a random system in the under- 
coupling regime cannot be approximated as a uniform 
slab with the average refractive index n e //, despite its 
quasimodes exhibit similar features as the Fabry-Perot 
modes formed by the reflections from the slab bound- 
aries. Since in our calculation the refractive index out- 
side the random system is set to n e ff, there would be 
no quasimodes if the random system were replaced by 
a dielectric slab of ??<•/■/'. Hence, the quasimodes in the 
under-coupling regime are not formed by the boundary 
reflection. In the over-coupling regime, the feedback from 
the scatterers deep inside the system becomes dominant, 
and the interference of multiply scattered waves lead to 
spatial localization of the quasimodes. 

Next we study the lasing modes in the random system 
with uniform gain and compare them to the quasimodes. 
rii is constant everywhere within the random system, so 
that the gain length l g — l/ki — l/k^rii in the dielectric 
layers is equal to that in the air gaps. Using the method 
described in the previous section, we find the frequency 
and threshold gain of each lasing mode. We calculate 
the lasing modes in the same random systems as in Fig. 
1 within the same wavelength range (500-750nm). The 
frequency k$ and normalized threshold ki/(ki) of each 
lasing mode are plotted in Fig. 1(c) for comparison with 
the quasimodes. It is clear that there exists one-to-one 
correspondence between the lasing modes and the quasi- 
modes for the random systems in both over-coupling and 
under-coupling regimes. For the system with rid = 2.0, 
the lasing modes match well the quasimodes, with only 
slight difference between ki/(ki) and koi/{koi) for the rel- 
atively leaky modes. For the system with rid — 1-05, 
the deviation of the lasing modes from the quasimodes is 
more evident, especially for those modes with large decay 
rates. Such deviation can be explained by the modifica- 
tion of transfer matrix M. In the passive system, fcoi is 
constant but ki = koin(x) varies spatially. With the in- 
troduction of uniform gain, ki becomes constant within 
the random system, and the feedback inside the random 
system is caused only by the contrast in the real part of 
the wavevector k r = kon(x) between the dielectric lay- 
ers and the air gaps. With a decrease in the scattering 
strength, k^ in the passive system gets larger, and the 
ratio of the feedback caused by the contrast in ki to that 
in k r increases. The addition of uniform gain results in 
a bigger change of M, as it removes the feedback due to 
the inhomogeneity of ki . Moreover, since there is no gain 
outside the random system, ki suddenly drops to zero at 
the system boundary. This discontinuity of ki generates 
additional feedback for the lasing modes. In the weakly 
scattering system, the threshold gain is high. The larger 
drop of ki at the system boundary makes the additional 
feedback stronger. To check its contribution to lasing, we 
replace the random system with a uniform slab of n e f f 
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while keeping the same gain profile. Since the real part of 
the refractive index or fc r is homogeneous throughout the 
entire space, the feedback comes only from the disconti- 
nuity of ki at the slab boundaries. We find the lasing 
threshold in the uniform slab is significantly higher than 
that in the random system with rid — 1.05. This result 
confirms that for the random systems in Fig. 1, the ad- 
ditional feedback caused by the fcj discontinuity at the 
system boundary is weaker than the feedback due to the 
inhomogeneity of k r inside the random system. However, 
if we further reduce rid or L, the threshold gain increases, 
and the feedback from the system boundary due to gain 
discontinuity eventually plays a dominant role in the for- 
mation of lasing modes. 

We also compute the intensity distribution I(x) of each 
lasing mode at the threshold. I(x) is normalized such 
that its integration across the random system is equal 
to 1. Such normalization facilitates the comparison of 
the lasing mode profile to the quasimode profile. In Fig. 
2(a) [(b)], I(x) of the lasing mode is plotted together with 
that of the corresponding quasimode. Although the las- 
ing mode profiles in Figs. 2(a) and (b) arc quite different, 
they are nearly identical to those of the quasimodes. For 
the localized mode in the random system with rid = 2.0, 
I(x) of the lasing mode does not exhibit any visible dif- 
ference from that of the quasimode in Fig. 2(a). For the 
extended mode in the system with rid = 1-05, the lasing 
mode profile deviates slightly from the quasimode profile, 
especially near the system boundaries. This deviation re- 
sults from the modification of the transfer matrix M by 
the introduction of uniform gain across the random sys- 
tem. The modification is bigger in the under-coupling 
system, leading to larger difference in the mode profile. 

Finally we investigate the lasing modes under local ex- 
citation. In particular, f{x) = 1 for \x — x c \ < L\/2, 
f(x) = cxp[— \x — x c \/L2\ for Li/2 < \x — x c \ < Li/2 + 
21/2) and f(x) = elsewhere. The lasing mode frequency 
fco, the threshold gain ki = fco^i, and the spatial pro- 
file I(x) are calculated with the method described in the 
previous section. I(x) is normalized in the same way as 
that of quasimode for comparison. As an example, we 
consider the same random structures as in Fig. 1 and 
introduce gain to the central region x c = L/2 of length 
L p = Li + 4L 2 = 8.84 + 3.03 = 11.87^m (marked by 
two vertical lines in Fig. 2). Figures 3(a) plots fco ar >d 
ki/(ki) for all the lasing modes within the wavelength 
range of 500-750nm. Comparing with Fig. 1, we find 
some quasimodes fail to lase under local pumping, no 
matter how high the pumping level is. The rest modes 
lase but their wavefunctions can be significantly modified 
by the particular local excitation. The two modes shown 
in Fig. 2 both lase under the local pumping configuration 
we consider. Their intensity distributions are plotted in 
Fig. 2. The mode in Fig. 2(a) is localized within the 
pumped region, and its spatial profile is barely modified 
by the local gain. In contrast, the mode in Fig. 2(b) is 
spatially extended and has less overlap with the central 
gain region. The intensity distribution of the lasing mode 
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FIG. 3: (Color online) (a) Normalized threshold gain ki/{ki) 
versus the frequency fco of lasing modes in the random systems 
with rid = 2.0 (solid squares) and rid = 1-05 (open circles) un- 
der local excitation (between the two vertical lines in Fig. 2). 
(b) Normalized frequency spacing Afco/{Afco) of neighboring 
lasing modes in the systems with rid = 2.0 (solid squares) and 
rid ~ 1-05 (open circles) under local excitation. 



differs notably from that of the quasimode. The expo- 
nential growth of I(x) towards the system boundaries is 
suppressed outside the gain region, while inside the gain 
region I(x) grows exponentially towards the ends of the 
gain region at a rate higher than that of the quasimode. 
These behaviors can be explained by the spatial variation 
of gain. Outside the pumped region, there is no optical 
amplification thus light intensity does not increase expo- 
nentially. Within the pumped region, the faster intensity 
growth results from the higher threshold gain for las- 
ing with local pumping than that with global pumping. 
Nevertheless, the close match in the number and spatial 
position of intensity maxima justifies the correspondence 
of the lasing mode to the quasimode. We repeat the 
calculation with many modes under the same pumping 
configuration, and find the weight of a mode within the 
gain region is often enhanced. To quantify such enhance- 
ment, we introduce a parameter d which is equal to the 
ratio of I(x) integrated over the pumped region to that 
over the entire random system. We compare the values 
of 8 for the lasing modes under local excitation to that of 
the corresponding quasimodes. For the mode Fig. 2(b), 
S is increased from 0.33 for the quasimode to 0.41 for the 
lasing mode, while for the mode in Fig. 2(a) d remains 
at 0.98. Thus the effect of local pumping is stronger for 
the modes in the weakly scattering system. This is be- 
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cause when scattering is weak the local gain required for 
lasing is high. The feedback within the pumped region 
is greatly enhanced, leading to the modification of mode 
profile. 

We also investigate the fluctuations in threshold gain 
and frequency spacing of lasing modes under local excita- 
tion. Figure 3(a) shows the lasing threshold fluctuation 
for the random system with rid — 1-05 is smaller than 
that with rid = 2.0. Since the number of lasing modes 
under local pumping is usually less than that of quasi- 
modes, the average mode spacing (Afcn) is increased. Fig- 
ure 3(b) plots the frequency spacing Afc of adjacent las- 
ing modes normalized to the average value (Afco). There 
is more fluctuation in the mode spacing for the random 
system with rid = 2.0 than that with rid — 1-05. Hence, 
with local gain the frequency spacing of lasing modes is 
more regular in the under-coupling regime than in the 
over-coupling regime. This result is similar to that with 
uniform gain. 

Although the local pumping enhances the feedback 
within the pumped region, the feedback outside the 
pumped region cannot be neglected. To demonstrate 
this, we calculate the lasing modes in the reduced sys- 
tems of length L p by replacing the random structures 
outside the gain region with a homogeneous medium of 
n e ff. The reduced system has uniform gain instead of 
the gain profile f(x) in the original system. The results 
are shown in Fig. 4(a) for the system with rid = 2.0 and 
in Fig. 4(b) for the system with rid = 1-05. The number 
of lasing modes in the reduced system is less than that 
in the original system under local pumping. In fact, the 
lasing modes are generally different, with only exception 
for a few modes localized within the gain region in the 
system with rid — 2.0. Moreover, the lasing threshold 
in the reduced system is higher than that in the original 
system with local gain. These differences are attributed 
to the feedbacks from the random structure outside the 
pumped region of the original system. It demonstrates 
the scatterers in the unpumped region also provides feed- 
back for lasing. By comparing Figs. 4(a) and (b), we find 
the difference in the lasing threshold between the origi- 
nal system under local pumping and the reduced system 
is smaller for the system with rid = 1-05 than that with 
rid = 2.0. It indicates the contribution from the scat- 
terers outside the gain region to lasing is reduced as the 
system moves further into the under- coupling regime. 

We note that local pumping introduces inhomogene- 
ity in the imaginary part of the refractive index, which 
generates additional feedback for lasing. To check its ef- 
fect, we simulate lasing in a homogeneous medium with 
the average refractive index n e //. The local gain pro- 
file f(x) remains the same. Only the spatial variation of 
ki(x) — kofiif(x) provides feedback for lasing. As shown 
in Figs. 4(a) and (b), the lasing thresholds are much 
higher than those in the random systems, even for the 
system with rid = 1.05. This result demonstrate that for 
the random systems in Figs. 3 and 4, the feedbacks for 
lasing under local pumping are predominately caused by 
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FIG. 4: (Color online): Threshold gain hi/{ki) of lasing modes 
in the (original) random system of length 24.1/im with lo- 
cal excitation in the central region of length 11.87/im (solid 
square), compared to the threshold gain of lasing modes in 
the reduced system of length 11.87/xm under uniform excita- 
tion (open circle) and the threshold gain of lasing modes in 
the homogeneous medium with n e f / under local excitation in 
the region of length 11.87/xm (a) rid = 2.0, n e ff = 1.3361, (b) 
n d = 1.05, n eff = 1.0168. 



the inhomogeneities in the real part of the refractive in- 
dex n(x) or the wavevector k r (x) = kon(x). However, a 
further reduction in rid or L p could make the feedback 
due to the inhomogeneity of fej(x) significant. 



IV. CONCLUSION 

We have developed a numerical method to calculate 
the quasimodes of ID passive random systems and the 
lasing modes at the threshold with either global or local 
pumping. We identify two regimes for the quasimodes: 
over-coupling regime [L > £) and under-coupling regime 
(L <C £)■ In the under-coupling regime the electric field 
of a quasimode grows exponentially towards the system 
boundaries, while in the over-coupling regime the field 
maxima are located inside the random system. The fre- 
quency spacing of adjacent modes are more regular in 
the under-coupling regime, and there is less fluctuation 
in the decay rate. The distinct characteristic of the quasi- 
modes in the two regimes result from the different mech- 
anisms of mode formation. In an over-coupling system, 
the quasimodes are formed mainly by the interference of 
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multiply scattered waves by the particles in the interior 
of the random system. In contrast, the feedbacks from 
the scatterers close to the system boundaries play a dom- 
inant role in the formation of quasimodes in an under- 
coupling system. The contributions from the scatterers in 
the interior of the random system to the mode formation 
are weak but not negligible. They induce small fluctua- 
tions in mode spacing and decay rate. As the scattering 
strength is increased, the feedbacks from those scatter- 
ers in the interior of the system get stronger, and the 
frequency spacing of the quasimodes becomes more ran- 
dom. 

In the presence of uniform gain across the random sys- 
tem, the lasing modes (at the threshold) have one-to- 
one correspondence with the quasimodes in both over- 
coupling and under- coupling systems. However, the las- 
ing modes may differ slightly from the corresponding 
quasimodes in frequency and spatial profile, especially 
in the under-coupling systems. This is because the in- 
troduction of uniform gain removes the feedback caused 
by spatial inhomogeneity of the imaginary part of the 
wavevector within the random system and creates addi- 
tional feedback by the discontinuity of the imaginary part 
of the wavevector at the system boundaries. As long as 
the scattering is not too weak, the quasimodes are only 
slightly modified by the introduction of uniform gain to 
a random system and they serve as the lasing modes. 
This conclusion is consistent with that drawn from the 
time-dependent calculations^^. Hence, the knowledge 
of the decay rates of the quasimodes, in conjunction with 
the gain spectrum, can predict the first lasing mode. Be- 
cause of the correspondence between the lasing modes 
and the quasimodes, the frequency spacing of adjacent 
lasing modes is more regular in the under-coupling sys- 
tems with smaller mode-to-mode variations in the lasing 
threshold. 

When optical gain is introduced to a local region of the 



random system, some quasimodes cannot lase no matter 
how high the gain is. The rest modes can lase but their 
spatial profiles may be significantly modified. Such modi- 
fications originate from strong enhancement of feedbacks 
from the scatterers within the pumped region. It in- 
creases the weight of a lasing mode within the gain region. 
Nevertheless, the feedbacks from the scatterers outside 
the pumped region are not negligible. Moreover, the spa- 
tial variation in the imaginary part of the refractive index 
generates additional feedback for lasing. As the pumped 
region becomes smaller, the number of lasing modes is 
reduced, and the frequency spacing of lasing modes is in- 
creased. In an under-coupling system, the regularity in 
the lasing mode spacing remains under local excitation. 
Our calculation results will help to interpret the latest 
experimental observation o 22 ' 23 of spectral periodicity of 
lasing peaks in weakly scattered random systems under 
local pumping. We note that the effect of local excitation 
can be significant in an over-coupling system if the size of 
the pumped region is much smaller than the spatial ex- 
tend of a localized mode or the spatial overlap between 
the pumped region and the localized mode is extremely 
small. Hence, caution must be exerted in using the decay 
rates of quasimodes to predict the lasing threshold or the 
number of lasing modes under local excitation. Finally 
we comment that the increase in the mode concentration 
in the gain region by local pumping have distinct physi- 
cal mechanism from the absorption-induced localization 
of lasing modes in the pumped region^. The former is 
based on selective enhancement of feedback within the 
gain region, while the latter on the suppression of the 
feedback outside the pumped region by reabsorption. 
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